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Abstract. Lapses of transmission phase in transport through quantum dots are 
ubiquitous already in the absence of interaction, in which case their precise location 
is determined by the signs and magnitudes of the tunnelling matrix elements. 
However, actual measurements for a quantum dot embedded in an Aharonov-Bohm 
interferometer show systematic sequences of phase lapses separated by Coulomb peaks 
- an issue that attracted much attention and generated controversy. Using a two- 
level quantum dot as an example we show that this phenomenon can be accounted 
for by the combined effect of asymmetric dot-lead couplings (left lead/right lead 
asymmetry as well as different level broadening for different levels) and interaction- 
induced "population switching" of the levels, rendering this behaviour generic. We 
construct and analyse a mean field scheme for an interacting quantum dot, and 
investigate the properties of the mean field solution, paying special attention to the 
character of its dependence (continuous vs. discontinuous) on the chemical potential 
or gate voltage. 
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1. Introduction 

Recent systematic studies [H El E] of current transmission through quantum dot (QD) 
embedded into an arm of an Aharonov-Bohm interferometer [H EJ El [7J , uncovered an 
unusual, correlated behaviour of transmission phase as a function of the gate voltage. 
Namely, between any two consecutive Coulomb blockade peaks the transmission phase 
suffers one abrupt change (phase lapse) of — tt. 

This surprising feature cannot be understood within the framework of a non- 
interacting QD model [3 [9l [TOl [11], where the presence or absence of a phase lapse 
between the two transmission peaks is determined by the relative signs of the tunnelling 
matrix elements coupling the corresponding QD levels to the leads. Roughly speaking, 
two adjacent peaks are separated by a phase lapse as long as the product of the four 
matrix elements, coupling each of the two levels to the two leads, is positive P, [TOT ITT]. 
Since experimentally there is no way to control these signs in a typical QD [12] . this 
would suggest an approximately 50% probability of the presence of a phase lapse between 
the two consecutive peaks, in disagreement with experimental data. This dictates that 
the Coulomb (charging) interaction between electrons in the QD must be accounted for 
at some level E3J El CEEl EE] . 

Earlier investigations [TST [TT] of the interaction effects in QDs resulted in the 
notion of "population switching" between the broad (strongly coupled to the leads) 
and narrow dot levels with varying gate voltage (or equivalently, with varying chemical 
potential). This phenomenon, which subsequently attracted attention of both theorists 
[HI HSl [201 12H [221 [231 12H [25] and experimentalists (see, e.g., reference [26]), consists 
in the narrow levels being shifted upward due to the Coulomb potential of electrons 
accumulated at a broad level, which remains near the Fermi energy ("hovers") over an 
extended range of the gate voltage/chemical potential values. Within the latter range, 
successive narrow levels from time to time get rapidly filled with electrons (thereby 
emptying the broad level, hence the term "population switching" [151 EH CHI H]) an d 
shift downward below the Fermi level. Available results [TFl ETJ [28] suggest that this 
switching may be either continuous or discontinuous, although no systematic study of 
the two scenarios has been performed. We note that already in references J27J [28] mean 
field approach (self-consistent Hartree approximation) has been employed. Apart from 
reference [17], where a QD with only one lead was considered, these early studies were 
all concerned with models where the absolute values of the tunnelling coupling of each 
individual level to the right and left lead were the same ( "left-right symmetry" ) . 

It is probably due to the latter circumstance that for a long time no attention 
has been paid to another generic interaction-induced mechanism which affects both 
the energy level structure and the transmission phase behaviour in an interacting QD. 
This mechanism, which is not effective in the opposite-sign left-right symmetric-coupling 
case only, in the context of more conventional solid state physics corresponds to forming 
excitonic correlations between different bands [29] . Within the mean- field approach to 
interacting QDs [30] . a similar scenario consists in forming the off-diagonal average 
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values between the QD states, leading to interaction-induced hybridisation between the 
dot energy levels. We note that the importance of considering the generic, right-left 
asymmetric dot-lead coupling values was pointed out independently in reference [20] . 

For the case of a strongly-interacting QD in the Coulomb blockade regime one finds 
that even a small right-left asymmetry typically results in a large effective hybridisation 
(again in analogy with conventional impurity problems [29]), leading in turn to a change 
of the coupling sign for the effective dot levels ("minus" changed to "plus"), and to a 
presence of phase lapse between the two transmission peaks. These findings, reported 
earlier in reference [30J, are in line with recent functional renormalisation group results 
obtained for both two-level [201 EH E2] and multi-level [21] interacting dots. We also 
note the recent perturbative treatment [23] and treatments based on the renormalisation 
group/Bethe ansatz approach [241 [25] . 

In the present article, we investigate transmission through a QD in the spinless 
case at zero temperature. We begin by describing the behaviour of transmission 
phase for multi-level non-interacting dot in section [21 The mean-field approach to 
an interacting two-level QD is formulated in section [31 The remainder of the paper 
is devoted to illustrating and discussing the two mechanisms, whose interplay brings 
about the abrupt changes of the transmission phase (as a function of the gate voltage or 
chemical potential) between the Coulomb blockade peaks. First, in section @] we discuss 
the "excitonic" restructuring of the QD spectrum, operational in the presence of a 
right-left asymmetry. Second, in section [5] we address the two scenarios (continuous and 
discontinuous) of "population switching" in the right-left symmetric case. The interplay 
between the two mechanisms is briefly described in section [61 The overall conclusions, 
along with some remarks on the generality of our results and a summary of outstanding 
questions, are relegated to section [71 

2. Phase Lapses in the Non-Interacting Case 

Transmission phase lapses in a multilevel dot in the absence of electron- electron interac- 
tion. Dependence on the dot-lead coupling strength. The role of dot-lead coupling signs. 
General expression for the transmission amplitude. Example: a four-level dot. 

In the present section, we discuss the behaviour of transmission phase in the case of a 
generic multi-level QD, in the absence of electron-electron interactions [9| [10 | fTT j [T3l [T^] . 
While our analysis here is far from being comprehensive, it allows to draw three 
important conclusions: 

(i) Phase lapses represent a generic property of transmission through QDs. In other 
words, the transmission amplitude t tr vanishing at certain values of the chemical 
potential or the gate voltage (at which point [TOl [14] the value of the transmission 
phase jumps by — 7r) does not impose any restrictive condition on the QD parameters. 
We quote reference [30] for an analysis of what is often perceived as a contradiction 
between this statement and the Friedel sum rule. 
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Figure 1. The model system, composed of a wire (chain) and a two-level dot, 
equations ([2H5]). Fully coupled dot corresponds to A = 1. The arrows show incoming 
and outgoing waves, cf. equation ([6]) 

(ii) As mentioned in the Introduction, in the non-interacting case the location of phase 
lapses with respect to transmission peaks is crucially affected by the relative signs of 
tunnelling matrix elements coupling the QD levels to the right and left leads [TUj 15]. In 
particular, a necessary (but possibly not always sufficient) condition for a single phase 
lapse to occur between the two successive dot levels (corresponding to the two successive 
transmission peaks) is that the product of the tunnelling matrix elements coupling the 
two corresponding levels to the two leads is negative (we stress that this can be severely 
modified in the presence of electron-electron interactions and an asymmetric QD-lead 
coupling, see sections 0] and [5] below). 

(iii) Beyond the above conclusion, the actual distribution of transmission zeroes with 
respect to the dot levels strongly depends on the values of the QD parameters, and this 
dependence can be rather complex. 

Other analyses of the transmission phase in the absence of interaction [8] have 
demonstrated that disorder or geometry alone would not give rise to sequences of 
correlated inter-resonance phase lapses. The latter may however result from a presence 
of a very broad dot level [3T] . 

The Hamiltonian of a generic non-interacting dot with Md levels E{ (see figure [1]) 
is given by 

M d 

n d = Y l {E i - l i)^d j . (i) 

i=l 

The dot is coupled to a one-dimensional wire 

= -9 £ ( a fe+i + a l+i a j) - a* £ d ) d v ( 2 ) 
z 3 3 

(where t is half bandwidth and the half-integer index j labels the sites in the wire) , by 
the perturbation V, 

V = Vt + g (c f i/2C-i/2 + cl 1/2 c 1/2 ) , (3) 



Phase lapses in quantum dots 



5 



Vn 



- 2 4 (4 } c_ 1/2 + ag } c 1/2 ) + H.c. , (4) 



where al (ag ) are the real tunnelling elements coupling the zth dot level to the left 
(right) lead; these are assumed to be small in comparison to the bandwidth in the wire, 
It. The second term on the r. h. s. in equation ([3D corresponds to cutting the link 
between sites j = — 1/2 and j = 1/2 of the wire. It is instructive to consider the system 
with varying perturbation strength, A, 

H x = H d + H w + XV (5) 

(figure [1]) . While for A = 1 this corresponds to a fully coupled dot (no direct hopping 
between right and left leads), characterised by a certain sequence of transmission 
peaks and transmission phase lapses, at A — > one recovers a featureless transmission 
amplitude, t tr (n) = 1 (hence Q tr = 0), of the decoupled case (unperturbed wire). As A 
increases from to 1, the profile of 0< r (/i) evolves between these two limiting cases. 

When searching for an eigenfunction ip which away from the dot has the form of a 
left-moving, partially reflected wave, 

!ti r exp(— ikx), x < — 1/2, 

(6) 
exp(— ikx) + run exp(ifczr), x>l/2, 

(where t tr and trr are transmission and right-right reflection amplitudes respectively) 
one has to solve a system of M^ + 2 linear equations for the quantities l/t tr , rjm/t tr , and 
the values ipi, i = 1, of the wave function on the dot levels (these correspond to 
the Schrodinger equation Tixip = —tip cos A; written for the wire sites j = ±1/2 and for 
the dot levels). Zeroes of the transmission coefficient t tr (at A ^ 0) obviously coincide 
with those of the principal determinant D of this linear system. It is straightforward to 
find that 



Docsinfc 2(1-A)t-A 2 E^^ IK e -^), (7) 




where e = — t cos k. We will now analyse this expression in the two limiting cases. 
A <C 1 - weakly coupled dot. As A decreases, the zeroes e = Zi approach the dot levels, 
e = Ei. To leading order in A, we find Z, = E i + \ 2 a ( £ ) aji /(2t). Thus, (i) the direction in 
which the transmission zero is shifted with respect to the corresponding dot level E^ is 
determined solely by the sign of the product of the tunnelling elements, <jj = signa^ ag\ 
We thus find that whenever all such signs are the same, each inter-level energy interval 
has exactly one transmission zero, whereas otherwise there arise intervals which contain 
either two such zeroes or none at all. (ii) Not counting the spurious zero at k = 0, there 
are exactly transmission zeroes. At A — > 0, each zero approaches the corresponding 
level Ei, cancelling its transmission peak and thereby restoring the featureless tt r (fj) 
[and tr (//)] at A = 0. 

A — > 1 - fully coupled dot. As the value of A increases, the transmission zeroes move 
further away from their corresponding dot levels. While at a finite A they obviously 
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cannot cross any of the dot levels to drift into neighbouring inter-level intervals, in 
those intervals which contain pairs of zeroes the zeroes may meet and disappear from 
the real energy axis; likewise, new pairs of zeroes may appear in some inter-level intervals. 
In addition, one or more zeroes which at smaller A may have been located below the 
lowest dot level or above the highest one may now move out of the conduction band. 
In the fully coupled case of A = 1, the r. h. s. of equation ([7]) is proportional to the 
polynomial of the power Md — 1, which guarantees that the maximal possible number 
of transmission zeroes is Md — 1; this number increases to Md for < A < 1 [it follows 
that at A — > 1 a transmission zero located outside the dot energy range must move out 
of the conduction band along the real energy axis]. When the two adjacent levels are 
characterised by the same signs <Xj = cr i+1 , it is certain that there is a transmission zero 
(or possibly an odd number of zeroes) between them; otherwise there may be either two 
zeroes (or, in principle, an even number of zeroes) or none. 

It appears that the possibility to have more than two zeroes in any inter-level 
interval requires a fine-tuning of the parameters of QD, and cannot be considered as 
generic. Barring these special cases, we see that had we compared this unrealistic non- 
interacting scenario to the experimentally observed correlated occurrence of exactly one 
transmission zero between every two adjacent dot levels, this would indeed require all 
the level coupling signs, o~i, to be the same. 



Furthermore, when at A = 1 for one of the dot levels Ei the quantity 



l (i) Mi 
K a R I 



is much smaller than the energy distances to the neighbouring levels, such a weakly 
coupled level is approached by a single transmission zero at 

-i 
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(8) 



We note that in this case, the direction from which the transmission zero would approach 
is not determined solely by o^. 

Finally, for reference purposes we quote the full expression for the transmission 
coefficient of a generic dot as given by equation (jSJ), 



tg>(fl,)=-iy/t*-f 



M d „(*)„(*) 



2(1-A)t + A 2 ^ 



o> a 
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i(2A - A - 2)tJt 2 -fi 2 - 



a ( L) 2 t + {afft - 2(1 - A)(/i - iy/W^J^a^a 



-a 2 E 
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At 



(9) 



Equation (Q contains the full information about the locations of phase lapses and 
transmission peaks. The latter are shifted with respect to the dot energy levels Ei, 
but the values of these shifts contain a pre-factor of the order of |a^o^ \/t 2 1 and in 
most cases of interest can be treated as small. 
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In particular, we find that at A = 1, and \o,^r\ ^ \Ei — E{\ for all i and /, the 
transmission phase Qt r (fi) in the vicinity of a dot level Ei is given by 

= ^ + r ^')' + '°gr ,n -do) 



(i) CO (0 (») x 2 

a R a L - a R a L 



V Y^l^i A{Ei-Ei)t 2 



This corresponds to the expected smooth increase of the transmission phase by tt with 
the value of /i increasing and sweeping E{. When the chemical potential lies away from 
the band edges, \fi\ < 1, the width of the step is of the order of the level broadening 



, l2 v t2 ~^> (11) 

and the sum in the denominator of equation ( fTOi) shifts the position of the centre of this 
step (coinciding with the transmission amplitude maximum) on the energy axis with 
respect to the bare dot level Ei. 

This complex behaviour of transmission phase is exemplified by figure [2] for the 
case of a four-level dot with o~i = {1, —1, —1, 1}. At a relatively small value of A = 0.3 
(dotted line), transmission phase remains close to zero except in the immediate vicinity 
of the dot levels, and there is a phase lapse of 7r near every level. With increasing A, 
the continuous increase of Qtr by rr in the vicinity of each level becomes progressively 
less steep, and the phase lapses shift further away from the levels. The directions 
and rates of these shifts reflect the differences in the coupling signs o~i and coupling 
magnitudes. While for all values of A there is exactly one phase lapse between E 2 and 
E 3 and none between E 3 and E±, the pair of phase lapses located between E\ and E 2 
approach each other, as shown by the dashed-dotted line (A = 0.75). They eventually 
merge and disappear, and with further increase of A a new pair of transmission zeroes 
emerges outside the dot energy range at \i < E\ (see the dashed line, A = 0.971). With 
increasing A, one of these new phase lapses moves towards E\, whereas the other moves 
rapidly to the left, disappearing in the fully coupled case of A = 1 (solid black line). 
The solid green line illustrates the effect of reducing the coupling of a single dot level E% 
to the leads in the A = 1 case. We see that the increase of transmission phase, Q tr (fj,), 

(3) 

by 7r near E 3 becomes steeper as we reduce a R L , and a phase lapse approaches E 3 from 

(3) 

the left, "annihilating" the smooth increase in the limit a RL — > 0. We also note that 
with decreasing coupling to the third QD level, the two of phase lapses located to the 
left of E\ and to the right of E± move away from the dot levels and disappear. 

In the general case of a non-interacting multi-level dot, the behaviour of 0t r (/x, A) 
reflects the analytical properties of the complex transmission amplitude t tr , which are 
probably not known in detail. In any case, these properties are far too cumbersome to 
try and analyse the effect of charging interaction on Q tr in such a generic case. One is 
therefore left with the option to consider the effects of interaction in a simple model case 

| This is given by 4I\ = ^[(ag) 2 + (a^) 2 ]p(l/2) 1 where p(l/2) is the local density of states at the 
terminal point of a lead, e.g., d(cj , 2 C\/2) / djj, at A = 1. 



Phase lapses in quantum dots 

E 




A, = 0.3 



X = 0.75 



X = 0.971 



A,= l 



X = 1, reduced a^ R 



Figure 2. (colour) The complex behaviour of the transmission phase for different 
regimes of QD-lead coupling. The dependence of &tr on /x for a non-interacting four- 
level dot with af/t = {0.05,0.07,0.1,0.08} and ct§ /t = {0.08,-0.08,-0.12,0.06}. 
The dot spectrum is given by E t /t = {0,0.03,0.07, .12}. The "perturbation strength" 
A (cf. figure [1]) is equal to 0.3, 0.75 0.971, and 1 for dotted, dashed-dotted, dashed, 
and solid black lines respectively. The green solid line corresponds to the fully coupled 
dot (A = 1) with reduced values of the tunnelling matrix elements for the 3rd level: 
4 3) = OMt and <4 3) = -0.0361 



in order to identify the relevant correlation-induced mechanisms with the hope that the 
results will prove generic at a qualitative level. 

3. The Mean Field Scheme for an Interacting Two-Level Dot 

The model for a two-level interacting QD. Rotation of the fermionic operators on the dot 
and the "intra- dot hopping" as a measure of right-left asymmetry. Mean field decoupling 
and mean field equations. Invariance of the mean field scheme with respect to the choice 
of the basis. Properties of the effective non-interacting system. 

In this section, we describe our method of calculation of the transmission phase, 
4r (jLt), for an interacting two-level QD. After introducing the model, we focus on the 
mean field decoupling and on the properties of the resultant non-interacting system. 

The minimal model for studying the phase lapse mechanism includes a two-level 

QD, 

H QD = (e{ 0) - n)d\dx + (4 0) - n)d\d 2 + Ud\d\d 2 d x . (12) 
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Here, the operators dj with % — 1, 2 annihilate electrons on the two dot sites (with bare 
energies {E^}, E 2 > E^). The QD is coupled to the two leads by the tunnelling 
term 

V T = ~d[ (a L C-i/ 2 + OflCi/ 2 ) - -d\ (b L C-i/ 2 + b R Ci/ 2 ) + h.c. (13) 

The operators dj (with half- integer j) are defined on the tight-binding sites of the left and 
right lead (cf. figured]). In the U = case, the location of the transmission phase lapse 
with respect to the dot levels is determined by the coupling sign, = signaROL&R&L; 
in particular, in the opposite-sign case of = — 1 the phase lapse occurs outside the 
inter-level energy interval. 

In spite of the simplicity of this model, (fT2HT3"|) . no exact analytical solution for 
general values of parameters has been found so far. In a recent paper [30], the present 
writers suggested a mean field approach to this problem. Here, we will further explore 
the generality of our mean field scheme and the physical properties of the mean field 
solutions. 

At the outset, we confine ourselves to studying the case when the values of 
tunnelling couplings in equation (TIB"]) obey the following condition: 

b 2 R -bl = al-a%, (14) 

which corresponds to a certain 3D subspace of the full space of all values of aL,R and 
b LR (the latter are assumed to be real). While we did not attempt an investigation 
of the case when the constraint ( |T4"j) is not satisfied, there is an expectation that no 
significant physics is missed by making this assumption (possibly barring a few isolated 
singular cases), which however offers an important technical benefit. Indeed, by making 
an orthogonal transformation of the QD fermion operators, 

(i 1 )=o(f), 0={ C °^ ~ sin M (15) 
\ d 2 J \ "2 / \ sirup cosf J 

with tan ip = (an — a^jKbi — b R ) we find that the tunnelling term, 
V T = - - {(az,cl 1/2 + a R c{ /2 ) (cos ipdx - sin ipd 2 ) + 

+ (bLcL 1/2 + b R c\ /2 )(sm(pd~i + cosy^)} + H.c. , 
then reduces to a simple form, 

V T = --a (cl 1/2 + c f 1/2 ) di - h) (c f _ 1/2 - c T 1/2 ) d 2 + H.c. (16) 

At the same time, the transformation ([To']) changes the form of the dot Hamiltonian, 
equation (fl2l) . to 

U QD = (E{ 0) - 11)^ + (Ei 0) - ix)d\d 2 - — {d\d 2 + dtdx) + Ud\d\d 2 d x . (17) 

Thus, the system can be formally viewed as a quantum dot with the "site energies" E^\ 
and the "intra-dot hopping" wq; coupling to the leads is now left-right symmetric, with 
the "coupling sign" a — — 1 (figure [3]). Further analysis will be carried out in terms of 
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Figure 3. The model system, composed of a wire (chain) and a two-level dot, 
equations (TIT))) and (fT?)) . 



this new QD, assuming £3 > -E 1 } • We note that the quantity 

1 (Ei 0) -E[ 0) )(b L -b R )(a L -a R ) 

Wo = — — — — (18) 

4 a L (a L - a R ) + b L (b L - b R ) 

is in reality a measure of left-right asymmetry of the original dot coupling, equation ( {TBI : 
by varying wq, one can probe both opposite sign (a^ = —1) and same-sign (a^ = 1) 
cases (see section H] below). 

The mean-field calculation entails decoupling the interaction term in equation ( fTTl) . 

d\d\d 2 di ^d\di(d 2 d 2 ) + d\d2 (^1^1) — {d[d\) (^2^2) — 

- dldzidldx) - d\d x {d\d 2 ) + I (d\d 2 )\ 2 , (19) 

which, when substituted in equation (fl7|) . yields an effective non-interacting dot, 

Hf F = (Ei - ^)<5di + (^2 - lAd\d 2 - ^(d\d 2 + dldi). (20) 

The self-consistency conditions take the form of three coupled mean-field equations, 

El = E{ 0) + U(d\d 2 ), E 2 = Ef ] + Uidl^), (21) 
w = w Q + 2U(d\d 2 ). (22) 

Here, the relevant average values are obtained from the thermodynamic potential, 
VLmfi of the effective non-interacting system with the Hamiltonian, 

KI = K F + n w + V T + t - (c{ /2 c_ 1/2 + cl 1/2 c 1/2 ) (23) 
[cf. equations (EH3])] according to 

(djdx) = d^Mp/dEr , (dld 2 ) = dn MF /dE 2 , (d\d 2 ) = -dn MF /dw.(2A) 

An exact and convenient method to evaluate Qmf is described in reference [30J. 

Clearly, the advantage of the choice (fT4"|) lies in the fact that at the mean-field level, 
interaction effects amount to a rather simple renormalisation of the coefficients in the 
single-particle part of the Hamiltonian, equation (flTl) . We note that a similar mean field 
scheme can be constructed for a symmetric same-sign QD, when the constraint (T141) is 
replaced by ai = a R , &l = b R . This case, which is analysed in reference [30J, will not be 
discussed here; we note also that a unitary transformation of the dot operators allows 
one to recast the corresponding Hamiltonian in the form (TX&HTTT) with b = 0. 
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It is important to notice that had we carried out the decoupling, equation (fT9"j) . 
in any other basis of the dot operators (including the original one, that of di )2 ), the 
resultant system of mean field equations would have been equivalent to ( 12X14221) . This 
means that the two operations: rotating the basis [cf. equation ({15]) ] and performing 
the Hartree - Fock decoupling [equation (|T9l) ] are commutative. Indeed, consider an 
arbitrary orthogonal transformation of the dot operators, 

(dx\_ Ql ( dj\ o , = /cos0 -sin0\ 



\ d 2 I \ d 2 / \ sin ip cos if> ) 

in the QD Hamiltonian, equation (fTTj) . Then in the new basis of operators d[ 2 we have: 

e[ 0) ' = cos 2 + Ef ] sin 2 - ^w sin 20 , (26) 

and similar equations for Ei and w f ; the interaction term remains unchanged. Now 
suppose that the "renormalised" quantities Ei j2 and w satisfy the mean field equations, 
f l2T]) . (1221) . and (12^1) . and do the same transformation O' after the decoupling (|T9|) . i.e., 
in the mean field Hamiltonian (1201) . yielding 

E[ = Ei cos 2 <p + E 2 sin 2 - -w sin 20 , (27) 



E' 2 = sin 2 + E 2 cos 2 + sin 20, (28) 



1 

2 

w' = (Ex - E 2 ) sin 2(f) + w cos 20 . (29) 

On the other hand, had we chosen to perform the decoupling in the d[ 2 basis, we would 
have obtained the mean field equations 

E[ = Ef ] ' + U dn MF /dE' 2 , E' 2 = E { 2 ] ' + Udn MF /dE[, (30) 
w' = w' - 2Udn MF /dw'. (31) 



It remains to verify that the quantities E[ 2 and w' as obtained from equations (12714291) 
solve the system ( |30ll3Tl) . Since 



8VLmf 9Qmf 2 j , 9Qmf . 2 , dVt MF . 

= — = = — cos 0H = — sm sin 20 (32) 

dE[ dEi 8E 2 dw Y y ' 

etc., this is easily done by a direct inspection. We conclude that the two systems of 
mean field equations, (12114221) and (130143T1) are equivalent to each other, hence the results 
are indeed independent on the choice of basis. This statement illustrates the fact that 
our mean- field approximation is a conserving one [32], and holds also if one replaces O' 
in equation (J25l) with an arbitrary unitary matrix. 

As we will see below, under certain conditions one finds that for a given value of 
/i, the mean field equations ([2114221) for Ei, E 2 , and w may have multiple solutions. In 
this case, one must choose the solution which corresponds to the lowest value of the full 
thermodynamic potential including the constant (non-operator) terms in equation (fT9|) . 

n = n MF - u{d\d 1 ) (d\d 2 ) + u{d\d 2 } 2 . (33) 
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While the invariance of Qmf under the transformation (I25p is obvious, it is 
straightforward to check that the sum of the last two terms on the r. h. s. also 
does not change, so that the entire quantity fl is independent on the choice of basis. 

Once the values of mean field parameters Ei >2 an d w are determined, the 
transmission phase Qtr [corresponding to the effective non-interacting model (1201) ] can 
be evaluated from 



2 



^tane„(,) = , + ^E 1 -,)+a\E 2 - : ) + 2,a%yt 

As a function of /x, Q tr shows two smooth steps of +n centred at the transmission peaks, 
/i = /ii j2 with 

/i 1>2 = (Eh. + E 2 )/2 T \[{E X - E 2 ) 2 + w 2 + Aa^/t 2 } 1 ' 2 . (35) 

More precisely, at \i = /ii j2 the quantity 6 tr /7r has half-integer values. Friedel sum rule 
then implies [30J that the same holds for the electron population change due to the dot, 
Ndot = N(/j,) — N w (fj), where iV(/z) is the total number of carriers in the system, and 
N w (n) is the number of carriers at an unperturbed (connected) wire alone [see equation 
(J2J)], evaluated at the same value of \x. 

The positions of these peaks /ii i2 are slightly shifted outwards with respect to the 
mean-field dot energy levels [eigenvalues of (1201)]: 

E lt2 = [E x + E 2 )/2 t \[{E X - E 2 f + w 2 } 1/2 . (36) 

Equation f|31|) determines Q tr up to a shift by a multiple of n. One can easily investigate 
the evolution of Q tr (/j,) for the effective non- interacting model (|20|) as the "interaction 
strength" A [cf. equation ([51)] varies from to 1. We thus find that this shift should be 
chosen in such a way that 

^ { -3vr/2 for a 2 < b 2 . (37) 

In addition, at \i = Z with 

E 2 a 2 - E x b 2 , s 

the transmission phase suffers a lapse of — tt [i.e., Q t r includes a term, —n8(fi — Z)\. 

In the range of parameters where multiple solutions to the mean field equations 
arise, it is possible that, e.g., \i = Z or \i = /i 12 corresponds to a thermodynamically 
unstable solution. This situation and its implications for the overall 0t r (/x) profile will 
be discussed in greater detail in sections [5] and El We shall now proceed with analysing 
the properties of mean field solutions for various values of parameters. 



4. First Mechanism for Abrupt Phase Change between Transmission 
Peaks: Effects of Left-Right Asymmetry And Excitonic Correlations 

Overview of mean field results: "phase diagram". Excitonic mechanism in "phase" 1: 
eigenstates of the dot are linear combinations of site states, hence effective coupling sign 
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-0.6 -0.2 0.2 ^ 0.6 1 

a 

Figure 4. (colour) The "phase diagram" of a two-level QD shown in figure [3] The 
parameters are U = O.lt, e[ 0) = 0, E { ° ] = 0.004*, and \/a 2 + b 2 = 0.125*. The axes 
represent the 1-2 level asymmetry, a = (|a| — |&|)/V» 2 + 6 2 , and the dimensionless 
intra-dot hopping, k = w q /[(e[ 0) - E { ® ] ) 2 + w 2 ] 1 / 2 . 



change due to off- diagonal correlations on the dot. Evolution of the effective dot param- 
eters with varying chemical potential \x. Excitonic mechanism becomes ineffective when 
approaching either the 1-2 or left-right symmetric situations. 

We begin our discussion of the mean-field results for a two-level QD with presenting 
a typical mean- field "phase diagram" (figure H]). This shows how the location of the 
transmission phase lapse Z with respect to the two transmission peaks [see equations 

and (|5S|) ] depends on the two dimensionless QD parameters, 

a= H - |&l , k = - (39) 

Of these, a is a measure of the 1-2 level broadening asymmetry, whereas k is related 
to left-right asymmetry of the levels coupling to the two leads. In this section, we 
will be interested in the generic situation when k is sufficiently large. In figure HJ this 
corresponds to "phases" 1 and 2 (blue and red). We will see that this typically gives 
rise to a hybridisation between the dot levels, which in turn results in the phase lapse of 
—7i occurring between the two peaks ("phase" 1). "Phase" 2, which occupies a narrow 
stripe near the 1-2 symmetric case, is characterised by the occurrence of a — 7r-phase 
lapse outside the region [i\ < fi < /i 2 . 

The physics underlying this correlated restructuring of the QD spectrum at larger 
k is thus somewhat similar to that of exciton formation in a semiconductor. A closer 
analogy can be drawn with the "excitonic correlations" between localised and extended 
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states in impurity problems • At smaller k — > we encounter another mechanism, 
that of "population switching" as discussed earlier for QDs in the Coulomb blockade 
regime [T5l [T7J ETJ [28] (section [5]). The entire "phase diagram", figure HI can be 
interpreted in the context of interpolation between these two regimes (section [6]) . 

Let us first formally examine the role of left-right asymmetry in our mean-field 
scheme, starting with equation (fT9]) . There, the diagonal and off-diagonal terms contain 
two types of quantum average values, 

"l,2 = (4,2^1,2) , "12 = (d\d 2 ) . (40) 

In terms of the "transformed" QD, equation (fl7|) . these correspond to the two dot site 
occupancies, and to the intra-dot "excitonic" [29J hybridisation respectively 

Let us first suppose that the value of wq is sufficiently small, so that in the absence 
of interaction the coupling of the two QD eigenstates (l36j) to the dot is opposite-sign, 
<7(°> = -l: 

< (2Z<°> - E{ 0) - E^ f - (E{ 0) - Ef ] f , Z®= 2 . (41) 

a z — cr 

In terms of figure HI this corresponds to the area below the thin solid line. We recall 
that in this situation, the non-interacting dot would have the phase lapse located outside 
the inter-level energy interval [although strictly speaking it may still barely fall between 
the two transmission peaks, which are slightly shifted with respect to the dot levels, 
equations (1551) and (1551) ]. 

As will be discussed in more detail below, the typical situation in the large-U case 
is that, due to large values of n X2 in equation (1221) . w is strongly renormalised in the 
relevant energy region around the phase lapse. If the value of w becomes sufficiently 
large, the effective dot sign o will change from —1 to 1, for this is obviously the case for 

\w\^\E 2 - E x \ (42) 

[when the product of the couplings of the two effective dot eigenstates, [d\ ± d 2 )/\/2, 
to the leads is given by (a 2 — 6 2 ) 2 /4 > 0]. The phase lapse will then be located between 
the two mean- field dot levels, E\ < Z < E 2 . 

This "excitonic" mechanism of the QD sign change is operational within the blue 
region of our "phase diagram" (figure HI "phase" 1), which is defined as the parameter 
region where the phase lapse of — tt occurs between the two transmission peaks and the 
properties of the effective dot vary continuously with //. We see that this blue region 
extends well below the line denoting the change of the original sign, o~(°\ which means 
that this situation is rather generic. It is further exemplified by figure [51 showing the 
evolution of the dot properties with \x for a particular choice of parameters. 

When the chemical potential lies well below the dot levels, the latter are unoccupied 
and the dot parameters are close to their bare values. In particular, the "intra-dot 
hopping" w [solid line in figure [5] (a)} is close to w = 0.0015£, which is small in 
comparison with the difference between the two "site energies" E^ 2 (dashed and dashed- 
dotted lines), E 2 — E\ w E^ — E^ = 0.004t. This ensures that coupling to the leads 
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Figure 5. Mean field properties of a QD with a = 0.3, k = 0.35, U = O.lt, E\ ' = 0, 
— 0.004i, and Va 2 + b 2 = 0.125i plotted vs. chemical potential \x. (a) "Intra- 
dot hopping" w (solid line) and the differences between mean field "site energies" 
and chemical potential, E\,2 — H (dashed and dashed-dotted lines). Dotted lines show 
the differences between the mean field energy levels and chemical potential, E\2 ~ M- 
(b) "Off-diagonal average" n.12 (lower solid line), "site occupancies" 711,2 (dashed and 
dashed-dotted lines), and the mean-field level occupancies ri\ t 2 (dotted lines). The 
upper solid line shows the transmission phase, Qtr (right scale). 



is opposite-sign, a = = — 1 [cf. equation (HIT) ]. The wave function of the lower QD 
eigenstate in this regime is mostly localised on the site 1 of the QD, which for our choice 
of a = 0.3 > is the one that is stronger coupled to the leads (a > b, see figure |3]); the 
lower level is therefore broader than the upper one [equation (1111) ]. 

With increasing /i, the population of this level [and hence the occupancy of site 1, 
dashed line in figure [5] (b )} begins to grow. The Coulomb repulsion term, U fii in the 
second equation (|2"T|) pushes the energy of the other site, E 2 , upwards, and its population 
fi2 [dashed-dotted line in figure [5] (b )} remains low. The energy E% (which is very close 
to the lower QD energy level, E\) eventually crosses the chemical potential, resulting in 
a smooth increase of fii and of the transmission phase tr [upper solid line in figure [5] 
(b)} in agreement with the Friedel sum rule [30] . 

While the average of the two dot energy levels, (Ei + E 2 )/2 [see equation (|36|) ] 
does not depend on w, the energy of the lower level Ei goes down when w increases, 
making such an increase energetically favourable in the partially-occupied regime of 
< "hi + fi2 < 2. We see that indeed the value of w begins to increase with increasing 
/i, which is accompanied by an increase of hybridisation [off-diagonal average value ^12, 
lower solid line in figure [5] (b)\. The nature of eigenstates begins to change continuously, 
and it is no longer possible to identify the lower eigenstate with site 1; at the same 
time, a large difference arises between the site energies Ei j2 and the energy eigenvalues 
£1,2 [dotted lines in figured] (a)]. While both the site energies Ei^ and site occupancies 
^1,2 eventually cross^jl, which might look reminiscent of the usual population switching 

§ We find fi\ — fi2 at [i ~ 0.08592;. The hybridisation reaches a maximum of fi\2 ~ 0.411 at fi ss 0.0852£, 
where fix — hi ~ 0.1. 
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scenario [T51 HZJ EZl ES] (see also section ED, the eigenvalues i?i j2 never coincide with 
each other, and the same holds for the respective level occupancies, 



[dotted lines in figure [5] (b)\. 

In particular, near the crossing point E\ = E 2 the level energies E\ >2 are pushed 
far apart by a large w, which peaks in this region, thereby satisfying condition fj42|) . 
Hence the dot becomes same-sign, a — 1, and when (also in this region of values of /x) 
the — 7r-phase lapse occurs with chemical potential crossing the transmission zero, the 
latter is located between the two level energies/transmission peaks. 

With a further increase of /i, the population of the QD increases towards fii+n 2 = 2, 
and the energy gain associated with large w becomes less pronounced. The value of w 
thus begins to decrease towards wq, and it is in this region that the second level crossing, 
E 2 = /i, occurs, accompanied by another smooth increase in ir . The site energies E lj2 
eventually cross again at fi ~ 0.145t, reverting to their original order. 

The presence of the = 1 area in figure H] (above the thin solid line) is indicative 
of the fact that our Hamiltonian, equation ( TT7T) . which is characterised by the opposite- 
sign "site coupling", a = —1, is suitable for probing both same-sign and opposite-sign 
original level coupling cases. We note that in figure HI most of the = 1 area is 
occupied by "phase" 1. 

The mechanism of interaction-induced coupling sign change becomes ineffective 
when approaching the line a = 0, corresponding to equal absolute values of coupling of 
the two QD sites to the leads. The reasons for this are two-fold: (i) a stronger increase 
in w/\E 2 — Ei\ is required to reach the coupling sign change in this case, (ii) in parallel 
with the usual "population switching" scenario, when the difference in the broadening 
of the two levels becomes smaller, larger U is needed to make the values of E\ and E 2 
cross. Thus if |a| is decreased while U is kept constant, the site energies E± j2 cease to 
cross, and in the partially-filled QD regime stay progressively further away from each 
other. This decreases the level overlap and hence the ability of the system to form n\ 2 
and thereby renormalise w. At the same time, the energy gain associated with larger w 
in the partially-filled region becomes smaller. 

Thus the coupling sign change does not occur for a relatively narrow "red" region 
( "phase" 2) around the a = line in figure HJ where the phase lapse of — 7r is located 
outside the energy interval between the two transmission peaks, fii < fi < fi 2 . "Phase" 
2 is also characterised by a continuous evolution of the dot properties with varying \i. 
The area occupied by "phase" 2 becomes smaller with increasing U or \wq\; it is located 
below the line denoting the sign change o^ -' = — 1 — > = 1. 

This "excitonic" mechanism also breaks down when approaching the u> = line. 
The reason for this is clear from figure [3j at w — 0, the contribution of the two virtual 
hopping paths (via the lead sites ±1/2) cancel each other, owing to the difference in the 
signs of coupling of site 2 to the right and left leads. Thus if we start with the wo = 
case, a non-zero off-diagonal average value ni 2 cannot be formed, and w remains equal 




(43) 
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to zero for all values of \i [see equation (I2"2"j) ]. Equivalently, in the w = case the dot 
sites 1 and 2 are coupled to even and odd combinations of electron wave functions in the 
two leads respectively; phases of these combinations remain fully independent of each 
other. 

It is, however, precisely in this region of small wo that the "population switching" 
mechanism becomes operational in its conventional form. In order to further understand 
the structure of our "phase diagram" , figure HI we now have to proceed with a more 
detailed analysis of the wo = case. 

5. Second Mechanism for Abrupt Phase Change between Transmission 
Peaks: Population Switching in the Symmetric Case — Discontinuous vs. 
Continuous Scenario 

"Population switching" in the right-left symmetric case. Evolution of the effective dot 
parameters in the discontinuous case. Multiple mean field solutions and the phase lapse 
renormalisation. Origins of the discontinuity. Continuous population switching in the 
absence of the Coulomb blockade. Effects of population switching and excitonic correla- 
tions on the transmission amplitude. Summary: excitonic correlations, continuous and 
discontinuous population switching. 

We will now consider the behaviour of an interacting two-level dot in the right-left 
symmetric case of Wq — 0. This situation was addressed earlier [15], [TTJ, [2TJ, [28| , and the 
associated notions of "hovering level" [15] or "population switching" [T7j were advanced 
in the literature. Nevertheless until recently (30J, [33] no clear distinction has been 
made between the two scenarios of continuous and discontinuous population switching. 
We will see that the difference between these two behaviours affects the magnitude 
of transmission phase lapse. More generally, the behaviour of a right-left symmetric 
QD turns out to be qualitatively different from the one found in the larger wq regime 
(coupling sign change due to "excitonic" correlations, section HJ). Once we clarify the 
effects of interaction at wo = 0, the entire "phase diagram", figure HI can be understood 
in terms of interpolation between these two regimes (section EJ). 

The mean-field analysis reported here has its obvious shortcomings. It is known 
( [20] |2T| |22| [231 EH EE], some of these references include approximate methods) that at 
least within a certain parameter range the discontinuity is smoothened (see also section 
[7] below). We include our mean- field results here as a convenient reference point for 
more elaborate analyses. We also note that available studies of the effects of quantum 
fluctuations [20] EH EH E3J EH [25] are not sufficient to conclude that the discontinuous 
evolution as described below is always an artifact of mean- field, to be " cured" in a more 
proper treatment. 

We already pointed out that in the wq = case the off-diagonal average value 
vanishes identically, hi 2 = 0, resulting in turn in the absence of effective intra-dot 
hopping, w = [equation ( 1221) ]. The "site energies" therefore coincide with the mean- 
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Figure 6. (colour) Mean-field properties of a two-level QD with right-left symmetry, 
plotted vs. chemical potential [i (black lines). QD parameter values are a = —0.6, 
w = 0, U = O.lt, E[ a) = 0, Ef ] = 0.004i, and vV + b 2 = 0.125*. (a) the differences 
between mean field energy levels and chemical potential, — fJ- (dashed and dashed- 
dotted lines), (b) Level occupancies, fix,2 (dashed and dashed-dotted lines). The solid 
line shows the transmission phase, Qtr (right scale). In both (a) and (b), black dotted 
and vertical solid lines denote the unstable solution and the discontinuous change of the 
QD state respectively, whereas the respective green lines correspond to the continuous- 
evolution case of smaller U and smaller |a|, viz. U = 0.03t, a = —0.25. 



field QD energy levels, E\ 2 = E X2 , and the same holds for site and level occupancies, 
^1,2 = n i,2 [see equations (|36j) and (j4"5j) ]. Since there is no hybridisation, the eigenstates 
of the QD do not change for all values of fi, and their respective couplings to the leads 
remain constant (no interaction- induced sign change can occur). Thus, the coupling of 
our QD to the leads remains opposite-sign, so that the transmission zero Z, equation 
( |38l) . always lies outside the energy interval between the mean field QD energy levels, 
on the side of the weaker-coupled level (i.e., for \a\ < \b\ we find Z < Ei^ if E% < E 2 , 
and Z > Ei t 2 in the opposite case of E\ > E 2 ). 

Typical behaviour of Expip) and 77, 12 (/i), as well as that of the transmission phase, 
tr (/x), is shown in figure El The figure corresponds to the a < case, when the coupling 
of the upper bare dot level (dot site 2) to the leads is stronger than that of the lower one, 
|6| > \a\ [equation (|39|) ]. We will now trace the changes of mean-field QD parameters 
with increasing fi, addressing first the case of stronger interaction effects (larger values 
of U and |a|), as shown in figure [6] by the black lines. We also refer to figure [7] for a 
schematic representation of the corresponding changes in the mutual orientations of the 
two dot levels £a,2, transmission zero Z and chemical potential /x. 
(i) When the chemical potential is well below the bare QD energy levels, fi < e[°1, 
the dot contribution to the overall density of states at the Fermi level comes mainly 
from the tail of the broader level E 2 . Therefore the occupancy of the dot site 2 [black 
dashed-dotted line in figure [6] (b)\ has a larger value, and with increasing /x increases 
at a higher rate, than that of site 1 (black dashed line). Owing to the large value of 
U, the mean field equations (!2~T|) then dictate that the mean field QD level energies Ei 
and E 2 [dashed and dashed-dotted lines in figure [6] (a)] eventually cross, E\ > E 2 for 
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Figure 7. (colour) Discontinuous population switching: schematic representation of 
the chemical potential dependence of the two level energies E 12 (dashed and dashed- 
dotted lines; it is assumed that the dot-lead couplings of levels 2 and 1 satisfy b > a). 
Dotted lines correspond to the unstable solution, whereas the discontinuity is marked 
by a bold vertical line. In the lower part, locations of -Ei,2, M an d the transmission 
zero Z are depicted schematically for five points A-E, showing explicitly the QD level 
rearrangement along with the switching in the location of Z. 
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/i > — 0.8t. The smooth increase of the dot population, and hence of the transmission 
phase Q t r{l^) [black solid line in figure [6] (b)}, continues beyond this crossing point, (see 
also figure [7J points A, B). 

(ii) . With further increase of fi, the lower mean-field level E 2 crosses the Fermi level at 
\i = 0.0052i, resulting in the usual feature (shoulder) in Q tr . While n 2 (fi) grows apace 
with Q tr , the value of na(/-t) remains low, passing through a maximum of ni w 0.017 at 
jj, — O.Oli. The transmission zero Z remains well above the Fermi level, Z > E\ > fi 
(figure d points B, C). At ji > 0.0075t, two additional solutions to the mean field 
equations appear [dotted lines in figures [6] (a) and (b)}. At the beginning, these new 
solutions are characterised by higher values of the total energy, fi(/x), which can be 
evaluated via equation (1331) . 

(iii) . As can be anticipated by considering the filled dot case 

with increasing /i the energy levels Ei j2 must eventually cross again, reverting to their 
original order (cf. figure [7J point E). In the present case, this second crossing occurs 
in a discontinuous manner, i.e., at some point it becomes energetically favourable to 
abruptly depopulate level 2 (which then shifts upwards) while putting most of the carrier 
population of the dot into level 1 (which is lowered). This is the "population switching" 
in its discontinuous form as proposed in reference [17] . Mathematically, with increasing 
\i within the multiple-solution region we eventually reach the point fi c ps 0.0275t, where 
the total energy values of the two lower-energy solutions cross. Hence a switch of the 
solution branch occurs, accompanied by a discontinuous change in all the QD properties 
(vertical solid lines in figure E]). While on approach to this critical value of fi from below 
the Fermi level was located below the transmission zero, the situation immediately 
following the transition is that of Z < E\ < fi < E 2 , i.e., the Fermi level is above the 
point of transmission phase lapse. The actual point Z = fi (and the associated phase 
lapse of —7t), however, is never crossed because it lies in the thermodynamically unstable 
region. This is illustrated by the dotted line showing transmission phases for unstable 
solutions in figure [6] (b), which includes a phase lapse of — tt at fi ~ 0.0381 The actual 
change of transmission phase at fi = fi c , 

AQ tr = -tt + nAN dot (44) 

includes both the contribution of this phase lapse and another term, related by the 
Friedel sum rule [30J to the jump of the "dot-related" particle number, Ndot = 
N(fi) — N w (fi). Here, N(fi) is the total number of particles in the system, whereas N w 
is the number of particles in an unperturbed (connected) wire, equation ([2]). Thus, Ndot 
includes both the dot occupancy, ni + n 2 , and the "dot-induced" change of population 
within the leads. N dot must always increase with decreasing gate voltage, or equivalently 
with increasing /x; in particular, it has a positive jump at \i — /i c , renormalising 
the value of phase jump A6 ir , equation (1441) . In our case, AQ tr ps —2.514, hence 
ANdot ~ 0.20 [note the difference of the latter from the jump in the dot level occupancy, 
A(ni +n 2 ) w 0.28]. In the schematic representation on figure the discontinuity (bold 
vertical line) occurs between the points C and D. 
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(iv) With further increasing /i, the (narrower) dot level 1 remains nearly filled, although 
its occupancy ni(/x) does not vary monotonously, passing through a local minimum of 
Hi ~ 0.987 at /x m 0.12t; the corresponding energy level E\ lies below the Fermi level. 
The occupancy of the other level, n 2 , increases with /i, with the level E 2 crossing the 
chemical potential at /x ~ 0.103t, resulting in another smooth increase of 0j r . We 
note that the same (broad) level E 2 crossing the Fermi level twice (above and below 
the jump) is a known feature ("hovering level") of the population switching scenario 
pj)J[T7]. The unstable solutions (dotted lines) provide a continuous connexion between 
the states of the system above and below the jump, with the level E\ crossing the Fermi 
level at /x ~ 0.089t (in the unstable region). With increasing /x, the pair of unstable 
solutions finally disappears at /i « 0.092t. In figured the second level crossing and the 
disappearance of the unstable solution occur between points D and E. 

The origins of discontinuity as found at /x = /x c become clear if one considers the 
case when one of the QD levels is fully decoupled from the rest of the system, e. g., 
a = (corresponding to a = —1) and Wq = 0. With increasing /x, the occupancy of 
the QD level 1 then changes abruptly form n\ = to ri\ = 1, at which point all other 
QD properties (n 2 , Ei t2 , 0t- r ) must suffer a jump as well. In the language of mean field 
(which is exact in this there are no quantum fluctuations of ni), this means the 

presence of multiple solutions in a finite region of values of /x near the jump point |17j . 
Indeed, assuming that the QD energy scales Ei 0) -E[°\ b 2 /t, and U are all much smaller 
than the bandwidth, 2t, we find that the value of £a(/x) is determined by a single mean 
field equation, 



where vq is the per-site density of states in the leads. It is easy to verify that there is 
at least one multiple solution region, where the system switches from the Ei > /x to 
Ei < ytx branch. These two branches of -&l(/x), connected by a segment of the singular 
line Ei(fi) = /x (where in the a = case the value of n\ is ill-defined) together form a 
z-shaped structure similar to that shown in figure [6] (or, in a cartoon form, in figures [7] 
and [9] below). When the value of a is increased from a = — 1 (corresponding to a ^ 0), 
this picture changes in a continuous fashion, so that in order to eliminate the multiple 
solution region (and hence the jump) altogether, a must exceed a certain finite value, 
a > a (U). In the opposite case of a < ao, the discontinuity persists, as exemplified by 
the a = —0.6, U = O.lt case shown in figure [6] (black lines) and discussed above. 

If in the latter case of a = —0.6, U = O.lt the value of a is increased further, one 
finds that the while the absolute value of the jump decreases (and the phase lapse value 
approaches — vr), the location of the discontinuity shifts further to the left. By the time 
the discontinuity disappears at a c {U) ~ .012, the phase lapse is outside the interval of 
the values of /x between the two transmission peaks. This situation changes when the 
value of U is smaller; we will now turn to the a = —0.25, U = 0.03t case, shown in 
figure [6] by the green lines. 

In this case, the mean field energy levels E\ and E 2 [green dashed and dashed-dotted 
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lines in figure [6] (a)} do not cross, E\ < E 2 for all values of /i. Nevertheless, they do 
come close to each other at fi ~ — 0.0025t, where E 2 — E\ « 0.0014t; a slight increase in 
| a | would give rise to a pair of crossings, E\ = E 2 > /i, in this region, without changing 
the overall picture. Throughout this low-/^ region (where the electron population on the 
dot, rii + n 2 , is below 0.5), the occupancy n 2 of the broader level 2 [green dashed-dotted 
line in figure [6] (b)} is larger than that of site 1 (green dashed line). The value of Ei 
crosses \i at \i ~ 0.007t; this is accompanied by a rapid increase in the value of ni, 
which exceeds n 2 for /i > 0.05t, giving rise to a sharp upturn in E 2 (fi) [cf. equation 
(121]) ]■ Around this point, it becomes favourable to depopulate level 2 (hence a downturn 
in n 2 ) while increasing n\. This process however happens continuously, exemplifying the 
scenario of continuous population switching, as encountered earlier in references [27] [28]. 
Comparing this scenario with that of U = O.lt, a = —0.6 discussed above, one may say 
that a slight non-monotonicity of n^ 2 (fi) noted in that case developed presently into the 
sharp maximum of n 2 at /i ~ 0.04t and absorbed the jump in the level occupancies. 

With a further increase of fi, level 1 remains nearly filled, whereas n 2 eventually 
starts to grow again, with E 2 crossing the Fermi level around fi ~ 0.031 This is 
accompanied by an increase in E\, with the difference E 2 — E\ approaching its bare 
value, E 2 — E^ = 0.004t, as the chemical potential increases beyond the dot energy 
range. 

Throughout the entire range of values of /j,, the transmission zero is located below 
the lower dot level, Z < E\. We note, however, that the points /i = /ii i2 as defined by 
equation ( [35]) (which we use as transmission peak locations, see below), are shifted with 
respect to those of the dot levels. This shift becomes relatively more pronounced when 
the dot levels lie close to each other (on the scale of the level widths), and if the value 
of a is not too large, the phase lapse may fall in between the two transmission "peaks" 
while remaining outside the interval between the two dot levels. This situation is realised 
in the present case, where the transmission phase [solid green line in figure [6] (b)} suffers 
a lapse of —it at fi ~ 0.005t (the point where n\ = n 2 ), shortly above the point of 
fi — fi\ ~ 0.002t. We note a strong asymmetry of the "peak", reflected in a rather 
irregular profile of 0£ r (/i) near the phase lapse. The rapid non-monotonous variation of 
the dot parameters, combined with a rather small inter- level distance E 2 — E\, results 
also in the absence of a well-defined "shoulder" of 0t r (/^) associated with the second 
transmission peak; the value of /i 2 is about 0.034t. 

These features are further illustrated by a plot of the transmission amplitude 
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as a function of /i (figured]), where the case of U = 0.03t, a = —0.25 is represented by 
the dashed line. We see that in this case, \t tr \ does not show the two well-separated peaks 




Figure 8. Transmission amplitude \t tr for an interacting two-level QD as a function 
of chemical potential: solid line, a — —0.6, U — O.lt, and wq = 0; dashed line, 
a = —0.25, U — 0.03i, and wo — 0; dashed-dotted line, a = 0.3, U — O.lt, and 
k = 0.35 (corresponding to wq = 0.0015t). The rest of parameters for all three cases 
are e[ 0) = 0, E^ 0) = 0.004i, and Va 2 + b 2 = 0.125*. The arrows show the locations of 
the corresponding transmission peaks as given by equation (|55|) . 



as anticipated in the Coulomb blockade regime. This is due to a relatively small value 
of U (and hence the small mean- field level separation). We also note that the profile of 
transmission is rather irregular. The points marked by the dashed arrows [locations of 
fj,x,2 as given by equation ( 1331) ] do not correspond to any particular features of the plot. 
This is not surprising, since these values of \i = n\p correspond to half-integer values 
of Ndot and should approach the transmission peak locations in the Coulomb blockade 
regime only. 

We note the presence of a transmission zero, corresponding to a phase lapse; this is 
surrounded by the two peaks reminiscent of "correlation induced resonances" which were 
first reported in reference [20J and subsequently explained by the Kondo-type physics 

[SUES]. 

The solid line in figure [8] corresponds to the discontinuous population switching 
case of U — O.lt and a = —0.6. We see the two well-separated Coulomb blockade peaks 
of roughly the same width, whose maxima lie very close to the values of equation 
( |38|) . As expected [301 there is no transmission zero, since the actual t tr {p) = 
point belongs to the unstable solution. The discontinuity is clearly seen on the plot. 
An interesting feature of both dashed and solid curves in figure [8] is the presence of two 
maxima where the value of |tt r (/i)| exactly equals one. An investigation of equation ( |46| ) 
suggests that this is always the case for the symmetric situation (w = 0), provided that 
the mean field value of \E 2 — E\\ (which in the partially filled regime is of the order of 
U) is larger than \2ab/t\ (here, we assume that both \E\ — E%\ and (a 2 + b 2 )/t are much 
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smaller than the bandwidth of the leads, 2t). 

This situation changes in the asymmetric, Wq 7^ exemplified by the dashed 

dotted line in figure El This line, which illustrates the effect of interaction-induced 
"excitonic correlations" (as discussed in section H] above) shows two clear Coulomb 
blockade peaks of unequal height, with a transmission zero in between. Again, equation 
( 155]) gives accurate values for the peak locations. 

* * * 

We are now in a position to summarise the behaviour of an interacting two-level 
dot in the three regimes considered thus far: 

(a) . Right-left asymmetry (w 7^ 0), excitonic correlations (section HJ). Within "phase" 
1 of figure HI this yields a phase lapse of —it between the two transmission peaks. 
The underlying mechanism is the interaction-induced change in the sign of the dot-lead 
coupling, the evolution of dot parameters with varying \i is continuous throughout, the 
transmission peaks are of unequal height and are separated by a transmission zero. 

(b) . Right-left symmetric (w = 0) case showing a discontinuous QD level population 
switching as a function of /i. The locus of the corresponding points in figure H] is at 
the lower edge of the green area marked as "phase" 3. While the interaction-induced 
sign change is impossible, the phase lapse is located between the two transmission peaks 
of height \ttr\ = 1, due to the re-ordering of respective locations of the peaks and the 
transmission zero at the point of discontinuity. The latter coincides with the phase lapse 
(whose magnitude is however reduced to a quantity between —it and 0) and with the 
transmission minimum (where the transmission amplitude retains a finite value). 

(c) Continuous population switching. This case would appear at a "phase diagram" for 
a lower value of U as the interval of the k = axis where the blue region of "phase 1" 
extends all the way down to this axis. While the example we considered lies outside the 
Coulomb blockade region, one observes a phase lapse of tt, coinciding with a transmission 
zero and surrounded by the "peaks" of transmission. This mechanism is driven by a 
strong non-monotonicity of the QD level occupancies in this regime. 

Scenarios (a) and (c) are both continuous and can be expected to be robust with 
respect to quantum fluctuations (not included in the present treatment). This may not 
be the case for the discontinuous scenario (b), which is therefore likely to be partly 
replaced with scenario (c). 

We now turn to the intermediate values of k in order to understand the structure 
of the entire "phase diagram" , figure HJ in terms of interpolation between these three 
regimes. 

6. Interplay between the Two Mechanisms: The Mean Field "Phase 
Diagram" 

Minimal value of 1-2 asymmetry required for the correlation-induced switching in the 
phase lapse location. Mean field "phase diagram" as a result of superposition of the two 



Phase lapses in quantum dots 



25 



correlated mechanisms: descriptions of different "phases". "Phase diagram" in the case 
of a smaller U and continuous population switching. 

In this section, we turn to the intermediate values of k [equation (1391) ] in order 
to understand the structure of the entire "phase diagram", figure HI in terms of 
interpolation between the different manifestations of interaction-induced correlations 
considered above - effective sign change due to excitonic correlations at larger k (section 
H]) vs. population switching at n = (section [5]). 

Both of these mechanisms result in the occurrence of a phase lapse between the 
two transmission peaks, and the underlying physics is partially similar in that once 
the chemical potential lies within a (broadened) dot level, the energetically preferable 
situation corresponds to the broader of the two QD levels being partially filled. This in 
turn implies that on approach of the value of fi to the dot level energies from the side of 
the weaker-coupled (bare) level (i.e., n < E x < E 2 ^ for a < 0), a level inversion must 
take place, with the broader level approaching the chemical potential first. At k = 
(no hybridisation) this inversion takes the form of an actual level crossing, whereas at 
k > only the "site energies" can cross, 

E^) = E 2 (fx c ) (47) 

(figure points A and B). Generally, the occurrence of the crossing, equation ff4T|) . 
indicates that the correlation effects are sufficiently strong to activate at least one of the 
mechanisms responsible for the occurrence of the phase lapse between the transmission 
peaks. In the case of chemical potential lying well below the dot energy levels, the two 
site occupancies are given by 

a 2 Jt 2 - a 2 _ b 2 Jt 2 - a 2 lA s 

n x « — ^ — , n 2 as — ^ *— 48 

irt 2 (E 1 - /i) irt 2 (E 2 -fj,) K ' 

[see reference [30], equations (16-17)]. In writing equation (148]) we assumed that 
t 3> Ei 2 — /U [which allows to use the constant value z/ = l/iji^t 2 — fx 2 ) for the density 
of states in the leaddjj and E X 2 — a 2 /t, b 2 /t. In addition, the value of w [equation 
(1221) ] should not be too large, w 2 < a 2 (t 2 - /i 2 )/t 2 and w 2 < b 2 (t 2 - /i 2 )/t 2 . Since in 
this range of values of chemical potential w stays close to its bare value, w ~ wo, [cf. 
figure [5] the latter is not a restrictive condition. The site energies crossing may occur 
below the transmission peaks only in the a < b (a < 0) case, which we will consider 
here. Then, equations (121]) yield the following condition for the crossing point fi c : 



UJt 2 -ii 2 c h 2 - a 2 
£ 1)2 ( Mc ) -fi e = goTI-W • ( 49 ) 

The two transmission peaks are well defined only if the dot levels are well separated 
at each peak. Hence the crossing point /i c must lie well below the lowest peak, in the 

|| The dependence of v$ on energy is indeed a weak effect, accounting for a slight asymmetry between 
a > and a < cases in figure SJ 
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region where the net occupancy of the dot, fix + n 2 , is small: 



2 T-,2 

- Mc > ri + r 2 = \A 2 - ^ > ( 5 °) 

where the quantity on the r. h. s. is the combined broadening of the two levels, cf. 
equation (fTTj) . We note that the condition fl50|) is also required for equations (jUD to 
hold. In terms of our 1-2 coupling asymmetry parameter a [see equation ( 139]) ]. equation 
( 1501) takes form 

7T f^ 0) - ^ 



[/ > t/ c « , = - y , (51) 



|a|V2^a 2 



where we wrote |a| instead of a in order to include also the a > b case, when the point 
\i c lies above the transmission peaks. The value of U c sets the scale for the interaction 
strength required to produce strong correlation effects. At a fixed U, these effects can 
be amplified by increasing the level coupling asymmetry so that 



1 - a 2 < 1 - oil 



7r»(Ef-Ef)* 



y — - ■ ( 52 ) 

For the values of parameters used in figure H] this yields a c ~ 0.09, which is a fairly 
accurate estimate for the width of the red area, occupied by "phase" 2. We already 
mentioned that in this "phase" the phase lapse occurs outside the energy interval 
between the two transmission peaks, due to insufficiently strong correlation effects. 

We will now turn to figure H] and summarise the mean-field properties of QD within 
the parameter region corresponding to each "phase". Typical 0t r (/x) profiles for each of 
the "phases" are shown in the main panel of figure [9j 

"Phase" 1 (blue). In this "phase", a phase lapse of —tt occurs between the two 
transmission peaks [marked by the two boxes on the ©t r (/x) plot]. There is no 
discontinuity in the evolution of other QD parameters with varying fi, as shown in 
figure M by a schematic /ii (/i) -/i vs. \i plot which does not have a multiple- solution 
region. As explained in section HI the underlying mechanism is that of an effective 
coupling sign change due to excitonic correlations, and the corresponding region is in 
the upper part of the "phase diagram" (larger k), away from the a = axis. 
"Phase" 2 (red). A phase lapse of —n occurs outside the interval between the two level 
crossings, fii^ifi) = H- The corresponding area is a narrow stripe around the a = axis, 
and at small or moderate values of k its width can be estimated with the help of Eq. 
( 1521) . Owing to an insufficient 1-2 level coupling asymmetry, correlation effects are not 
strong enough to cause either a sign change or a discontinuity (note the absence of a 
multiple solution region in a schematic plot of fi\(fi) —jj)- Away from the a = axis the 
"site occupancies" (//) may still show a non- monotonous behaviour and continuous 
population switching (with the phase lapse located away from the transmission peaks). 
"Phase" 3 (ereen). In this "renormalised" phase lapse, AQ tr > —it, occurs 

between the two transmission peaks. This "phase" occupies the area adjacent to the 
k = axis, excluding the vicinity of the fully symmetric case, k = a = 0. The 
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Figure 9. (colour) Typical behaviour of Ot r (/i) in different "phases" ( left; plots 
shifted vertically for convenience). Relative positions of transmission peaks, fii t 2, 
(boxes; also in the left panel) and the — 7r-phase lapses (circles) in "phases" 1-6 are 
clarified by the schematic plots of /xi — /z around the multiple-solution region (absent 
for "phases" 1-2). Solid (dashed) lines correspond to stable (unstable) solutions. The 
abrupt "switching" of solutions (vertical solid line) may either renormalise the phase 
lapse (when the — 7r-lapse lies in the unstable region) or result in a positive jump of 

e tr . 



underlying mechanism (discontinuous population switching) was discussed in section 
|5j Note the presence of a "fold" on the schematic plot of — fi, indicating the 

presence of a multiple solution region. The unstable part of the solution is shown by 
the dashed line and includes the transmission zero (circle), whereas the vertical solid 
line corresponds to a discontinuous jump between the two stable branches. Since at 
k = the multiple solution region has a finite width, the effects of right-left asymmetry 
(favouring the change of the coupling sign) cannot eliminate the discontinuity also at 
small but nonzero values of k, corresponding to this "phase" . 

The case of a fully decoupled level a = ±1, k = is essentially the one considered 
in the reference p2] for a multi- level QD. We note that the value of phase lapse at these 
two points also differs from — n, contrary to earlier expectations [17]. For the values 
of the QD parameters used in figure HI we find AQ tr ~ —2.4 at k = 0, a = — 1, and 
AQ tr ~ -2.5 at k = 0, a = 1. 

"Phase" 4 (cyan). This "phase" lies below the bold line which in figure H] forms the 
lower boundary of the "phase" 1 area and corresponds to the onset of discontinuity 
in the evolution of the dot parameters with /i. With decreasing k, this onset occurs 
in a continuous fashion, via the point where at some value of no the derivative of the 
QD parameters [such as dni(n)/dn] becomes infinite. With a further decrease of k this 
derivative changes sign [in the case of dfii(fi)/dfi, becomes positive], and a multiple 
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solutions area forms around Initially, this area is narrow and the unstable solution 
includes neither transmission zero nor transmission peaks (since in general fiQ dose not 
correspond to either of those). This is represented by the corresponding schematic 
plot of — \i in figure [H the jump between the two branches yields a discontinuous 
increase of the transmission phase, A0^ > 0, which is clearly visible in the 0t r (/i) plotj^jl 
to the right of the transmission zero, AO^ = — n. When the value of k is decreased 
further, the unstable part of solution spreads to include the transmission zero, and we 
find ourselves within the "phase" 3. Since "phase" 4 interpolates between "phases" 1 
and 3, it is clear that both phase jumps lie in the region between the two transmission 
peaks. 

"Phase" 5 (yellow). This "phase" occupies the areas adjacent to those of "phase" 3 
from the side of smaller values of the 1-2 level coupling asymmetry, |a|. It lies in the 
discontinuous region of the "phase diagram" below the bold solid line. The decrease of 
| at | within "phase" 3 causes the discontinuity (which in this case includes a jump "over" 
the transmission zero) to shift towards one of the transmission peaks. Eventually the 
unstable area moves to include this peak as well, as can be seen from the schematic plot 
of — \i. Thus the transmission peak is circumvented by a jump, as reflected by the 
absence of the corresponding box (and the corresponding inflexion point) on the Q tr (fjj) 
plot, which also includes a single renormalised phase lapse, AQ tr > —n. 
"Phase" 6 (magenta). When the value of \a\ in "phase" 5 is decreased further, the 
multiple solution area moves further away from the centre of the energy interval between 
the two transmission peaks. Eventually the unstable part (which still includes the 
transmission zero) clears this interval altogether, restoring the transmission peak to the 
stable branch. This corresponds to crossing the boundary from "phase" 5 into "phase" 
6. The renormalised phase lapse, AQ tr > —n, is then located outside the area between 
the two peaks. 

"Phase" 7 (not shown due to the small area it occupies). With a further decrease of 
| a J , the multiple solution area shrinks (while shifting away from the transmission peaks), 
and the transmission zero returns to a stable branch. The result is similar to "phase" 
4 above, with the only difference that the two phase jumps now lie outside the region 
between the two peaks. 

A further decrease in \a\ results in the disappearance of the multiple solution region 
and crossing the bold line into "phase" 2 so that the sequence of "phases" 5, 6, and 
7 interpolates between "phases" 3 and 2. The seven "phases" discussed above do not 
exhaust all the possibilities of mutual overlaps between the multiple solution region of 
values of \x (and associated discontinuity), transmission zeroes, and transmission peaks. 
Some of the other "phases" occupy minute areas (not shown) near the crossing of the 
boundary between "phases" 1 and 2 and the bold line; still others do not arise for the 
values of parameters used in figure HI 

We see that the whole of "phase diagram" , figure H] can be understood in terms 

\ When plotting Qtr(^) in this case, we used different values of the QD parameters in order to make 
this behaviour more pronounced. 
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of interpolation between the three cases: that of "phases" 2 (where the interaction 
does not affect the location of the phase lapse between the transmission peaks), 3 
(discontinuous population switching), and 1 (effective coupling sign change). Of these, 
the latter occupies the largest area. We will discuss the implications of this findings in 
the following section. 

It should be emphasised that in figure H] the boundaries between "phases" 1 and 
2, and between "phases" 3, 5, and 6 are a matter of convention and do not correspond 
to sharp transitions of any kind. Rather, they merely mark the changes in mutual 
locations of transmission peaks as defined by equation K35\) on one hand, and phase 
lapse/discontinuity on the other. 

When the value of U is reduced to U < 0.04* (while keeping the other QD parameter 
values in figure H] const ant), the bold line marking the onset of the discontinuity no longer 
intersects the 1-2 "phase" boundary. In this case, "phases" 5, 6, and 7 diasppear and 
the area of "phase" 1 formally extends down to the k = axis at a certain range of 
values of a. In the U = 0.03* case [see the green lines in figure [6], this range is given 
by ai ] « -0.31 < a < a ( f ] « -0.2 and c4 +) w 0.18 < a < a{ +) « 0.3. At this values 
of a and for small right-left asymmetry k < 1, the "continuous population switching" 
scenario as discussed in section formally results in an occurrence of the phase lapse 
between the two transmission peaks [as defined by equation (135]) ]. We stress however 
that at least for the moderate values of k, the entire region, a[ ^ < a < the QD 
is outside the Coulomb blockade regime (transmission peaks are not well separated, as 
exemplified by the dashed line in figure EJ). As a consequence, equation ( 135]) (used by us 
to define "phase" 1) looses accuracy, making our convention for distinguishing between 
"phases" 1 and 2 problematic. 

We note that the values of a[ , marking the onset of the discontinuous behaviour, 
are in a good agreement with equation (152]) . which at U — 0.3* yields \a c \ ~ 0.3. With 
a further increase of \a\ towards \a\ = 1 the Coulomb blockade effects set in even at 
U = 0.3*. 

7. Discussion 

Conclusions of the present mean-field study. Anticipated results for a mean-field treat- 
ment of multi-level dots. Role of quantum fluctuations: available results and outstanding 
questions. 

Broadly speaking, our mean-field treatment of the interacting two-level dot yields 
the following result, which holds provided that the interaction is sufficiently strong, see 
equation ( T5TT) : 

Irrespective of the original sign of the dot-lead coupling, the transmission phase lapse in 
an interacting QD generally occurs between the two transmission peaks. 

There are two distinct correlation-induced mechanisms which bring about this 
uniform situation. Of these, one is related to the off-diagonal ( "excitonic" ) correlations 
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on the dot [30] and requires the presence of a finite left-right asymmetry in the original 
dot-lead coupling (as can be expected generally in the experimental realisations). The 
other, which in its pure form is operational in the left-right symmetric case, has to do 
with the correlation- induced "population switching" [T5l [T7j . which can occur either 
discontinuously (in the Coulomb blockade regime) or continuously (this leads to the 
Coulomb blockade being lifted). 

At the most basic level, the two mechanisms share the same origin, familiar from 
the standard solid state physics: namely, when a band or an impurity level resides at the 
chemical potential and is therefore partially filled, energy may be gained by increasing 
the width of the band or by broadening the impurity level. In case of the "excitonic" 
mechanism, this "broadening" of the narrower QD level is achieved via a strong increase 
of hybrydisation with the broader level. This in turn is associated with the change of 
the coupling sign, causing the — 7r phase lapse to occur between the two transmission 
peaks. In case of the discontinuous population switching [17], the levels of a partially 
filled dot are abruptly rearranged in such a way that the narrower level never actually 
crosses the chemical potential. This "jump" of the narrow level (say, E\) from E\ > fi 
to E\ < fi with increasing fi is accompanied by a similar jump of the transmission 
zero Z and therefore leads to an abrupt change of a transmission phase, A<d tr > —n 
(renormalised phase lapse [30]). Continuous population switching in the relevant regime 
(which arises only for smaller values of U, see sections |5] and El) involves the two levels 
approaching the chemical potential at the same time (along with the transmission zero) 
and is accompanied by a phase lapse of —it in the absence of the two well-defined 
transmission peaks. 

The behaviour of a QD for general values of the left-right and 1-2 coupling 
asymmetries can be understood in terms of superposition of these two mechanisms. 
The continuous population switching scenario can evolve into the "excitonic" one via a 
smooth crossover with increasing left-right asymmetry. The interplay of these two with 
the discontinuous population switching mechanism, on the other hand, gives rise to a 
number of different intermediate "phases" . These are characterised by different numbers 
(1 or 2) and mutual locations of phase jump(s) and transmission peak(s). Our results 
suggest that on the whole, the behaviour found in most cases is the one dominated 
by the "excitonic" mechanism, with a phase lapse of — n located between the two well 
separated Coulomb blockade peaks. 

From the theoretical standpoint, our results give rise to the following two questions: 
(i) how are these findings generalised in the case of a multi-level interacting QD? (ii) 
What is the role of fluctuations, neglected in our mean-field treatment? We will now 
address these issues in some detail. 

(i) Multi-level dots within the mean-field approach. In general, one can expect that 
our conclusions will hold for a multilevel dot, with the two levels nearest the chemical 
potential playing the role of an "effective" two level QD, which at least at the mean 
field level would behave in a qualitative agreement with our results. One expected 
change is that the parameter area corresponding to the effective sign change ( "excitonic" ) 
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mechanism will be expanded at the expense of the other "phase diagram" regions 
(dominated in the mean field approach by the population switching). Indeed, for a two- 
level dot the ineffectiveness of this mechanism in the absence of the right-left asymmetry 
(section H]) can also be viewed as originating from the fact that each of the dot levels 
is coupled to a different "subset" of carriers in the lead (odd and even wavefunctions), 
making the inter-level hybridisation impossible. For a multi-level dot, the number of dot 
levels exceeds the number of such subsets (which of course remains equal to two, even 
and odd). Hence even if the two levels adjacent to the chemical potential have different 
coupling signs and only a small right-left asymmetry, each of these will hybridise with 
other levels further away (some hybridised combinations may even be decoupled from 
the leads [31] ■ This will in turn give rise to a left-right asymmetry of the "effective" 
two-level QD provided these other levels have asymmetric couplings to the two leads. 
On the other case, we saw (section H|) that even a moderate amount of initial left-right 
asymmetry is sufficient to activate the "excitonic" mechanism of the effective coupling 
sign change. 

These conclusions are in a qualitative agreement with the numerical results for 
multilevel dots [21]. We note, however, that these calculations [21] were performed 
using the functional renormalisation group method, and therefore include fluctuations 
at some level. 

We also note that within the mean- field approach the discontinuous evolution of 
the QD parameters (driven by the discontinuous population switching mechanism) will 
persist within the corresponding range of parameter values also in the case of a multi-level 
QD. The physical reason for this is the same as in the two-level case (section [5]) : suppose 
that one of the QD levels is fully uncoupled from the leads (and from the other levels). 
It is then clear that with increasing /i it will eventually be filled in a discontinuous 
manner (the occupancy jumping from to 1), leading to a discontinuity in all the 
QD properties p2]. Within the mean-field treatment, such a disontinuity originates 
from a jump between different solutions to the mean field equations, and therefore 
indicates a presence of a multiple-solution area in the parameter space. This area has a 
finite size, and therefore removing it altogether and thus eliminating the discontinuity 
requires a finite (as opposed to infinitesimal) dot-level coupling. Like in the two-level 
case, the discontinuous evolution will be accompanied by a renormalisation |30J of the 
corresponding phase lapse values, A9 ir > — n; the conductance at the corresponding 
values of chemical potential or gate voltage would not vanish. As described in section 
|6l there will also arise a number of borderline "phases" ; in some of these, the additional 
positive phase jumps would appear. 

(ii) Role of quantum fluctuations (validity of the mean field approach). The present mean 
field treatment does not include the effects of fluctuations. It is therefore important to 
understand to what extent do these alter the overall picture. While considerable effort 
has been made recently in this direction [20l [211 1221 [23J, [HJ [25], certain questions still 
remain unanswered. Below we will attempt to summarise the available results while 
pointing out the open problems. 
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First, we note that the available numerical and renormalisation group studies 
suggest that the generic effects of electron- electron interaction on the dot include the 
appearance of the phase lapse between the two transmission peaks. This is in qualitative 
agreement with the mean field results. On the other hand, the essentially many-body 
features like the "resonances" surrounding the transmission zero [201 El] , which originate 
from the Kondo-type physics [211 [25] cannot be reproduced at the mean field level. 

Our mean field results suggest the presence of two distinct correlated mechanisms 
causing the phase lapse to occur within the inter-level energy interval. This is in line with 
recent renormalisation group results [25], suggesting that the QD behaves differently in 
the right-left symmetric ("parallel effective field" [25J) and asymmetric ("tilted effective 
field") cases. 

Of the two mechanisms identified within the present mean field approach, the 
"excitonic" one (involving off-diagonal correlations on the dot, sectionH]), which requires 
the presence of some right-left asymmetry in the dot-lead coupling (w ^ 0), is found to 
be more generic. Indeed, it gives rise to the prevalent "phase" 1 in our "phase diagram" , 
figure HI This mechanism is not related to instabilities of any kind, and can be expected 
to remain robust beyond the mean field. While this is again in line with the available 
results, a more qualitative comparison can and should be made in order to confirm that 
we correctly identified the underlying physics. To this end, one should verify that the 
off-diagonal average, riyi = (g?!^) indeed shows a sharp peak in the regime of partial 
QD occupancy [cf. figure [5] (b)}. While the quantity n 12 should be readily available 
from numerical calculations, we are not aware of any published results for io We note 
that the diagonal average values, nip, have been calculated recently by the numerical 
renormalisation group |2l] and functional renormalisation group [22J methods, as well as 
analytically [25] . In the relevant range of values of parameters, their dependence on the 
gate voltage shows strong non-monotonicity (first noticed in the mean-field studies of 
the right-left symmetric case [271 EE]) and looks rather similar to the mean field results 
as seen in figure [5] (b ). Another interesting question is related to the degree of 1-2 level 
coupling asymmetry required for theis mechanism to be effective. The corresponding 
mean-field result, equation (152]) . differs from the one obtained via the renormalisation 
group calculation in the "local moment" regime [reference [25J, equation (70)] and a 
systematic numerical investigation is required in order to establish whether (and when) 
either of these is close to the actual value. 

The other mean field mechanism which gives rise to a phase lapse between the 
Coulomb blockade peaks is that of discontinuous population switching (section [5]). While 
the physical origins of the discontinuity are quite clear (see above), the associated 
notion of multiple solution region is restricted to the mean field approach. Thus, the 
discontinuous scenario will surely be strongly affected by the fluctuations. At present, 
it is not clear whether it survives in the exact treatment. Neither the discontinuity, nor 

+ Our results are in a qualitative agreement with Equation (52) of reference [25] which suggests that 
77.12 reaches a maximum in the general area of the point where h\ — fi2- This result |25j should become 
exact in the large- U, hi + fi2 — > 1 limit ("local moment regime" [25]). 
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the associated "renormalised" phase lapses [30] were found (except at isolated points 
|22j) in either the functional renormalisation group J2TJ [22], numerical renormalisation 
group [2TJ[22], or analytical renormalisation group/Bethe ansatz [2H [25] treatments. 

Lastly, we wish to turn to the situation of perfect right-left symmetry (w = 0) 
and argue that it represents a singular case where the effects of fluctuations are most 
pronounced. Indeed, owing to the symmetry of this case the hybridisation, n 12 , and 
hence the effective intra-dot hopping w remain equal to zero for all values of chemical 
potential or gate voltage (sectionHJ). Hence the transmission zero, equation ( |38i) . remains 
outside the energy interval between the two transmission peaks at all times (at least in 
the Coulomb blockade regime, which is of interest to us here). We conclude that within 
the mean field approach for w = the only possibility to observe a phase lapse between 
the two well separated (Coulomb blockade) transmission peaks is via a discontinuous 
restructuring of the spectrum, as explained in section [5] (figure [7j). In this case, the 
phase lapse is renormalised [30] , and is not accompanied by a transmission zero [301 S3] • 

On the other hand, the numerical results[2T], as well as those of the renormalisation 
group approach [Ml 125]. reliably indicate that at least for some values of 1-2 asymmetry 
in the right-left symmetric case, the phase lapse of exactly — n, associated with a 
transmission zero, does occur between the two Coulomb blockade peaks. This cannot be 
reconciled with the mean field picture at all, via any sort of renormalisation of the mean 
field parameters (which could be held responsible for the suppression of discontinuous 
behaviour at w ^ 0). Rather, fluctuations in this case must be giving rise to a 
totally new physics (albeit perhaps only in the narrow interval of values of chemical 
potential/gate voltage around the transmission zaro). We note that from the point 
of view of renormalisation group analysis [25] the wq = ("parallel effective field") 
situation does correspond to a special case, and its relationship to behaviour at w ^ 
remains unclear. 
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